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Abstract — Being inspired by pliasor analysis in linear circuit 
theory, and its algebraic counterpart - the AC-(operational)- 
calculus for sinusoids developed by W. Marten and W. Mathis - 
we define a complex structure on several spaces of real-valued 
elementary functions. This is used to algebraize inhomogeneous 
linear ordinary differential equations with inhomogenities stem- 
ming from these spaces. Thus we deduce an effective method to 
calculate particular solutions of these ODEs in a purely algebraic 
way. 
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I. Introduction 

When considering linear time-invariant systems, i.e. (at least 
in this paper) when considering linear ordinary differential 
equations with constant coefficients, a number of ad-hoc- 
methods are used to calculate the response to certain kind 
of input functions. In this regard, phasor analysis, the AC- 
calculus for sinusoids developed by W. Marten and W. Mathis, 
and "Ansatze" according to the special form of the input come 
to mind. 

We will show that all of these methods fit under a common 
heading - a heading, which we call the "full AC-calculus". 
This calculus can be characterized on one hand by its tendency 
to group certain functions into linear spaces of functions. Thus 
we are able to use the full force of linear algebra. Another 
characteristic of this calculus is the extensive use of complex 
structures, and the quite astonishing fact that a number of 
spaces of real-valued functions allow the introduction of a 
multiplication with complex scalars, so that they carry the 
additional structure of a complex linear space. This quality 
allows for nice formulae and elegant algorithms. 

This paper is still tentative, with proofs having been left 
out. Furthermore we do not yet show aU possibiUties of the 
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full AC-calculus. Nevertheless the results are sound, and we 
hope to give our readers a glimpse at what might be possible. 

II. Notation and Review of known facts 

Let / be the interval [0, oo). The set of n-times continuously 
differentiable real-valued functions on the interval / will be 
denoted by (7" (/) . We remind the reader of the fact that C" (/) 
is an infinite dimensional real linear space, i.e. superposition 
holds in C"(/). 

Given x G C'"'{I) we define a linear operator L on C"{I) 
by setting 

L{x) ~ anX^"'> + an-ix^"^'^^ H \- aix + a^x (1) 

with Oo, . . . ; a„ G R. We call the operator L normalized if 
a„ = 1. To L we assign the characteristic polynomial pL, 
which is 

Pl(s) = a„s" + a„_is""^ H \- ais + uq. (2) 

A inhomogeneous linear (ordinary) differential equation can 
now be written very conveniently as 

L{x) = r, (3) 

where r is a continous function on /. The corresponding 
homogeneous linear ODE is given by 

L{x) = 0. (4) 

Furthermore for x G C'^{I) and to G we define the vector 

x(to) as 

x(to):-(a;(<o),i(io),...,2:("-i)(to))^. (5) 

A typical initial-value problem can be written in the form 

L{x) = r, x(0) = xo (6) 

where r is like above and xo G R". 

The following theorem (cp. [1], §20, Theorem 1, §19 I. and 
V.) summarizes the known facts about the solutions of (|4]i, ^ 
and I©. 

Theorem 1: The solutions of (|4|i form an n— dimensional 
R-linear subspace X of If A is a real zero of the char- 

acteristic polynomial with multiplicity k, then the functions 

e^*,te^*,...,i'=-^e^* (7) 
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form a fc — dimensional subspace of X. If A + juj,X,uj G 
R, is a complex zero of the characteristic polynomial with 
multiplicity k - thus, because the coefficients oq, . . . , a„ are 
real, A — juj is a complex zero of multiplicity k, too - then 
the functions 

e^* sin(wi), ie^* sin(wt), . . . , t'^'^e^* sin(wi), (8) 
e^* cos(wt), ie^* cos(wt), . . . , i'^^^e^* cos{ujt) (9) 

form a 2/c — dimensional subspace of X. Consequently, 
the zeroes of the characteristic polynomial pl counted with 
multiplicity determine n R-linear independent solutions of 

The solutions of the inhomogeneous linear ODE (O are 
given by 

x = Xp+Xh (10) 

where Xp IS 3, fixed solution of (O and Xh is any solution of 
the homogeneous linear ODE (|4]i. Thus the solutions of (O 
form an affine space A, which is given by 

A^Xp+X. (11) 

Finally the initial-value problem (|6]l uniquely determines a 
"point" xq G a, i.e. a function a;o G C" {I), which then can 
be written in the form 

xo = Xp + Xho. (12) 
with a suitable function x^q ^ X. ■ 
In the language of system theory we refer to r as the input, 
to the operator L as a linear time-invariant system and to the 
solution Xq of the initial-value problem dSjl as the output or 
response. 

Now we can formulate the time-honoured general strategy 
to solve the initial- value problem ©I 

1) Find a particular solution Xp to the inhomogeneous linear 
ODE ©. 

2) Determine the zeroes of the characteristic polynomial pL 
with multiplicities, and calculate the R-linear indepen- 
dant solutions Xhi, ■ ■ ■ , Xhn to (|4]i according to (|7ll, (EJ 
and (ig. 

3) Make an "Ansatz" 

n 

X := Xp + '^^hiXhi, (13) 

i=l 

determine the derivatives and solve 

the system of inhomogeneous linear equations resulting 
from setting x(0) = xq with respect to the bi. 

Steps 2) and 3), while sometimes necessitating annoying 
and not very easy to do calculations, conceptually do not pose 
any problems. The main difficulty of the above algorithm is 
the task set up in 1): finding at least one particular solution to 
a given inhomogeneous linear ODE. 

Indeed in theoretical electrical engineering, there are several 
very well known approaches to this problem. All of them stem 
from the fact that, although the solution to the initial-value 
problem (|6]l is a unique function, its decomposition into a 
sum of two functions as in ( fT2] i is by far not unique (cp. [2], 
Chap. 4.3, in particular the introductory paragraph of section 
4.3.2). We will substantiate this remark by two examples. 



A. Zero-State and Zero-Input Response ([2], Chap. 6.1) 
In the situation given by ^ the unique solution of 

L{xm)^0, xm(0)-xo (14) 

is called the zero-input solution/response to the initial-value 
problem, the unique solution of 

Lixpo)^r, xpo(0)-0 (15) 

is called the zero-state (solution/response). Clearly, x^o is a 
solution of the homogeneous ODE. By the Unearity of L we 
have 

L{xpo + Xho) = L{xpQ) + L[xm) = r + = r, (16) 

and 

[XpO + Xm](0) = Xpo(O) + Xm(0) = + Xo = Xq. (17) 

Thus XpQ + x/io is the solution to (|6]). 

The advantage of this particular decomposition is seen, 
when the Laplace-transform is brought into play and thus we 
are able to define the powerful tool of the network function of 
an LTI-cii-cuit (cf. [3], Chap. 10, sect. 4.4). 

B. Transient and Steady-State 

A decidedly different decomposition is used, when we 
demand that the right-hand side r of ^ for some G R 
satisfies 

r G C^, (18) 

where €^ C C"(/) is defined by 

■■= { f : I ^ 'R-.t ^ acos(wi) +/3sin(wi) la,/? G R}, 

(19) 

and all the zeroes of the characteristic polynomial pl{^) 
have negative (< 0) real parts. Under these assumptions, 
there is a bounded function Xpr, which is a solution to the 
inhomogeneous ODE 

L{xpr) = r (20) 

and which is uniquely determined already by (l20b . Further- 
more we have 

Xpr G ^Ci; ■ 

(21) 

The function X IS called the steady-state solution of The 
corresponding homogeneous solution Xhr, which satisfies 

Xq — Xpr ^~ Xfir (22) 

is referred to as the transient solution. 

The method to determine X /p/p IS known as sinusoidal steady- 
state or phasor analysis (cf. [3], Chap. 9). While this method 
is considered by some as theoretically unsound (cf. e.g. [4], 
Chap. 3.3., where it is called a "rule of thumb"), the work 
of Marten and Mathis [5], [6], [7], [8], [9], as will be ai-gued 
below, has shown that this is not the case. 
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C. Steady-State response vs unbounded input 

While the above cited resuh on the existence and uniqueness 
of the steady-state solution can be generalized to inputs r G 
Cb{I), i.e. to arbitrary bounded continuous input functions 
[9], it has become customary in Control Theory to consider 
unbounded input and speak of the resulting response, which 
usually is unbounded, too, as a "steady-state", as well (cf. [4], 
Chap. 3.1 and 3.4). We point out the fact, that this linguistic 
lapse can be corrected [10], but we will not do so in this paper. 
Even worse, from this point onward, we will make the same 
use of the word "steady-state" and will call any solution Xp 
of the equation L{xp) = r a steady-state-solution, keeping in 
mind that in general, while existence is guaranteed for any 
r G C{I) by theorem [T] uniqueness usually is not. 

D. Some more notation 

We close this section by fixing the notation for some more 
sets of input-functions. Let A,ijJ G R, m G N U {0}. With 
R[i], we denote the R-linear space of polynomial functions. 
We define 

<tx:^ [t^ ae^*- : a^H}, (23) 

^x+juj -.^ {t^ acos{ujt)e^* + /3 sin(wi)e^* : a,/3 G R}, 

(24) 

gqj^ :={t^ p{t) ■ e^* : p{t) G R^ }, (25) 

e^'^ -.^{t^ p{t) • e^* : p{t) G R[t],deg{p{t)) < m }, 

(26) 

{ t p{t) ■ cos(wi) + q{t) ■ sin(a)i) : 

p{t),q{t)eR[t]}, 
(27) 

C:<P" := { i p{t) ■ cos(wi) -I- q{t) ■ sin(wt) : 

p{t),q{t) G R[t],m&x{deg{p{t)),deg{q{t))} < m}, 

(28) 

'S'Pa+j^ ■■={t^ Pit) cos(w<)e^' + q{t) sin(wt)e^* : 

p{t),q{t)eR[t]}, 
(29) 

and finally 

'^^x'+juj ■■={t^ p(t) cos(wt)e^* + q{t) sin(cjt)e^* : 

p{t),q{t) G R[t],max{deg(p(t)),deg(g(t))} < m}. 

(30) 

It is obvious that all of the above sets form R-linear spaces 
of real-valued functions - C^P^^ and C^Pa+ji^ being infinite 
dimensional - and we clearly have the inclusions 

e:^ c g:C' c c'p^, (3i) 

€a+,. C C e:*P,+,^ (32) 

and 

C^P™ = (B<Jii]l, (34) 

- e^^^. (35) 



For A = ti) = we have, that is equal to the space R[t] 
of polynomial functions. 

In addition, we have that each of these spaces is closed 
under differentiation, thus L induces a linear operator on each 
of them. The following sections are dedicated to a closer study 
of the action of L in each of the above cases. 

III. Complex Structures on the above function 

SPACES 

Earlier we have stressed the fact, that all of the above 
sets of functions are real linear spaces. Marten and Mathis 
have shown in the eighties that carries the structure of 
a complex linear space, as well, and from this were able to 
build up the AC-calculus for functions in C^^. In this and the 
next sections, step-by-step we will generalize this result to 
the space (S^PA+jw '^hus setting the foundation for the "full" 
AC-calculus. We will further look at how differentiation and, 
consequently, the action of the operator L fit together with this 
complex structure. 

A. AC-Calculus for Functions in vs. Phasor Analysis 

In this section we closely follow the line of thought set up 
in [5], [6], [7], [8] and [9]. 

Let x{t) — acos{Lot) + /3sin((jji) G 't^ and a + jb G C. 
Then we can define a multiplication with complex scalars 
on by setting 

j Q x{t) = j Q {a cos{u!t) -\- P sin{ujt)) 
:= —asm{ujt) + /3cos(w<) 

and continuing this rule by 

(a + jb) x{t) 
= {a + jb) Q {acos{Ljt) -\- (3sm{ujt)) 
:= {aa + bf3) ■ cos{ujt) -\- {a(3 — ba) ■ am{Ljt). (36) 

With this multiplication at hand, (i^ becomes a 1-dimensional 
complex linear space of functions. If we fix as a basis the 
function 1 1-^ cos{ujt), then we have 

a cos{ujt) + fi sin(ti>t) = (a — jfi) cos(wt). (37) 

Let us point out the fact that, if we consider a sinusoid 

t ^ z{t) ^ Acos{ujt + ^) (38) 

with A, a;, $ G R, then due to the addition rule for the cosine, 
we have 

Aco&{ijjt + $) = A - (cos($) cos(wt) - sin($) sin(wi)) 
= [A ■ (cos($) + j sin(<I>))] cos(wt) 
= (Ae^'*) 0cos(wt). (39) 

Thus the phasor associated to the sinusoid z(t) is the com- 
plex scalar with which the basis-function cos{u>t) has to be 
multiplied according to to get z{t). 
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We now look at the effects of differentiation on functions 
in €uj- For x{t) as above we find 



dt 



x{t) 



{a cos(wt) + /3 sin((jj(t))) 



dt 

= Lijj3 coB{Ljjt) — Ljasm{Lo{t)) 

= (ju) {a cos{Lut) + /3sin(w(i))) (40) 

= Uto)Qx{t). (41) 

Thus differentiating a function x(t) G Ci^; is the same as 
multiplying x{t) with the complex scalar juj according to 0. 

Now as an easy consequence, we have for aU functions 

X G ^(jj 

L{x) = plUuj) Q X (42) 

where pl again denotes the characteristic polynomial of the 
operator L. Thus we get 

Theorem 2: Let r G Cuj- If jtu is not a zero of the 
characteristic polynomial pL of the differential operator L, 
then the function 

1 

0r 



PlUuj) 

is a particular solution to the inhomogeneous linear ODE 

L{x) — r. 



Thus - if there had been any doubts - phasor analysis now 
is completely rehabilitated. 



B. A slight generalization 

We are now looking at the space <B\+j^. A typical function 
X in this set is given by 



x{t) — acos{Lot)e +Psm{ujt)e 



(43) 



with a, /3 e R. 

Here again, we can endow Ca+ji^ with a complex structure 
by setting 

{a + jb)Qxit) 
= (a + jb) (acos(wt)e^' + /?sin(tjt)e^*) 
:= {aa + &/3) cos(u;t)e^* + (a/3 - ba)sin{Lut)e^\(44) 

Ca+jo) thus becomes a 1-dimensional complex space, as well. 
The differential operator ^ acts on 2;A+ja; via 



dt 
d 

di 



x{t) 



{a cos(wt)e^* + (3 sm{ujt)e^^) 



At\ 



= {aX + Poj) cos(wt)e^* + {PX - auj) sm{ujt)e^\45) 

= (A + juj) {a cos{Lut)e^* + (3 sm{ujt)e^^) 

= {X + juj)Qx{t). (46) 

With the same reasoning as in the above case, we thus attain 
the following result: 



Theorem 3: Let r G ^\+juj- If A + ju) is not a zero of 
the characteristic polynomial pL of the differential operator 
L, then the function 

1 



r 



PL{X+juj) 

is a particular solution to the inhomogeneous linear ODE 

L{x) = r. 



IV. The full AC-Calculus 

The fuU AC-Calculus is a generalization to both phasor 
analysis, as described above, and the classic "Ansatz according 
to the right hand side". Thus it rehes heavily on Unear algebra. 

A. Polynomial input and input from 

Let us fix TO G N. We know, that is an m + 1- 

dimensional real linear space of functions. For the duration, 
we take as a fixed basis in €^P™ the functions {uk ■ fc G 
{ 0, . . . , TO }}, where we set 



Uk '■ t 1-^ t 



(47) 



As we know from theorem [T] those functions are linearly 
independent and, indeed, form a basis. 

2;*P™ is closed under differentiation, thus -4 gives a linear 



operator on 2;*^™ 
basis functions 
d 

d-t""' 



which is determined by its action on the 

d 



dt 
k-t 

kuh 



fc-i 



e 

1 + Xuk 



for fc > 1 and 



dt 



Uq 



dt 



Xe 



At 



Xt'' ■ e^* 



Xuq. 



(48) 



(49) 



If we identify the function 



x{t)^yakt''e^' 



k=0 



fe=0 



with the vector 



x= (ao,ai,...,Q;„)^, (50) 

then the action of ^ on given by matrix multiplication 

as 



/A 



dt 



\ 



\ 



TO 



X (51) 



where empty entries denote 0. From this it follows, that the 
action of the operator L is given by an (m+1) x (TO+l)-matrix 
as well. The next theorem can be shown by induction. 
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Theorem 4: Let again denote the characteristic polyno- 
mial of L. Then 

L{x) ~ iaki)o<k,i<m X 
where the matrix {aki)o<k,i<m is given by 

for k > I, 



PL (A) for k = I, 

[Qpt'\X) forfc<L 



(52) 



A slight exercise in matrix computation shows that this matrix 
{aki) is invertible, iff A is not a zero of pl- Thus we have 

Theorem 5: If r £ C;^™ is given by 



(53) 



fe=0 



and if A is not a zero of the characteristic polynomial pL, then 
a particular solution x to the equation L{x) ~ r is given by 



(54) 



where 

and successively 
1 / 



PL (A) 



•^m — l 



rm-l 



E 

k—7n — l-\-l 



k 

m — I 



(56) 

for 1 < I < ra. ■ 

Setting A = in the above discussion settles the case of 
polynomial input. 

B. Input from (L^'^ 

First we show how to make £^P™ into a complex linear 
space. For 

x{t) = p{t) cos{ut) + q{t) sin(cjt) (57) 

with real polynomials q{t) G R[t] and a, 6 G R we set 

(a + jb) x{t) 
= (a + j5) {p{t) cos{Lut) + q{t) sin(wi)) 
:— {a ■ p{t) + b ■ q{t)) cos{ujt) 

+ {a-q{t)-b-p{t))smiujt). (58) 

While is a 2(m + l)-dimensional real linear space, 

with this new scalar multiplication it becomes an m + 1- 
dimensional complex linear space. Fixing as a C-basis the set 
of functions {vk ■ < k < m}, where 



Vk 1 1-^ t*"' cos{u)t) 



(59) 



we see that differentiation acts as a C-linear operator on a; G 
£?P™. The corresponding matrix is given by the formula 



dt 



joj 2 

juj 3 



v 



juj m — 1 

ju m 



X. (60) 



In analogy with the above discussion we get 
Theorem 6: If r G £?P™ is given by 



"(i) = t'^jak cos{u!t) + Pk sinjut)) 

fc=0 
m 

= ^ 7fc t'' cos(wi) 



k=0 



with 7fc — ak — jPk for < fc < m, and if juj is not a zero 
of the characteristic polynomial pL, then a particular solution 
x to the equation L{x) = r is given by 



c{t) =^XkQt'' cos{ujt), 



fe=0 



(55) where 



(61) 



(62) 



and successively 
for 1 < I < m. 



E 



A;— m — ^+1 



- / 



{juj)xk 



V. Conclusion 

Obviously the full AC-Calculus, as demonstrated above, 
extends to the space SJ^A+jw' too, since we can define a 
multiplication with complex scalars for these functions, as 
well. 

Furthermore in this paper, we have not yet regarded the 
question, what will happen, when we have resonance, i.e. when 
in theorem |5] the scalar A or in theorem |6] the scalar joj is a 
multiple root of the characteristic polynomial p^. Keen readers 
will have guessed that the spaces 2;^a ^^x+jw ^re 

perfectly suited for handling these cases. In fact, the formulas 
given in theorems |5] and |6] can easily be adapted. We will leave 
this task for a later time. 

For now, we hope that we were able to convince the 
gentle reader that it is perfectly alright and algorithmically 
advantageous to consider sets of real functions as complex 
linear spaces. 
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Note added to the Electronic Version 

In this electronic document, some small typographical errors of 
the printed version were corrected. This especially refers to the last 
formula given in Theorem |6] 

Furthermore, for the convenience of the reader an abstract, key- 
words, MSG classification, and a short GV according to IEEE 
standards have been added to the arXiv-version. (Sept. 17th, 2007) 



